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Abstract 

We present the results of new Suzaku observations of the Coma Cluster, the X-ray brightest, 
nearby, merging system hosting a well studied, typical giant radio halo. It has been previously 
shown that, on the western side of the cluster, the radio brightness shows a much steeper gra¬ 
dient compared to other azimuths. XMM-Newton and Planck revealed a shock front along the 
southern half of the region associated with this steep radio gradient, suggesting that the radio 
emission is enhanced by particle acceleration associated with the shock passage. Suzaku 
demonstrates for the first time that this shock front extends northwards, tracing the entire 
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length of the western edge of the Coma radio halo. The shock is detected both in the temper¬ 
ature and X-ray surface brightness distributions and has a Mach number of around M 1.5. 
The locations of the surface brightness edges align well with the edge of the radio emission, 
while the obtained temperature profiles seem to suggest shocks located 125-185 kpc further 
out in radius. In addition, the shock strengths derived from the temperature and density jumps 
are in agreement when using extraction regions parallel to the radio halo edge, but become in¬ 
consistent with each other when derived from radial profiles centred on the Coma Cluster core. 
It is likely that, beyond mere projection effects, the geometry of the shock is more complex than 
a front with a single, uniform Mach number and an approximately spherically symmetric shape. 

Keywords: galaxies: clusters: individual: Coma — X-rays: galaxies: clusters 


1 Introduction 

Clusters of galaxies grow through mergers with other (sub-)elusters and by aeereting material from the 
surrounding large-seale strueture filaments. Major mergers between systems of eomparable masses 
release tremendous amounts of energy into the intraeluster medium (ICM). These events drive shoeks 
and turbulenee in the ICM, thereby heating the thermal diffuse gas, aeeelerating partieles to relativistie 
energies, and amplifying the intra-eluster magnetie fields (see reviews by e.g. Markevitch & Vikhlinin 
2007 ; Briiggen et al. 2012 and referenees therein). 

Evidenee for these mergers ean be readily identified at a variety of wavelengths. In the hot, 
diffuse ICM, merger shocks create discontinuities in the temperature, density and pressure profiles 
that can be measured through X-ray observations. More recently, shocks have also been detected using 
the Sunyaev Zel’dovich (SZ) effect, whose strength is directly proportional to the thermal pressure of 
the X-ray emitting plasma (e.g. Planck Collaboration et al. 2013 ). 

In addition, particles accelerated to relativistic energies during the mergers emit synchrotron 
radiation in the radio band. Depending mainly on their morphology and location with respect to the 
X-ray peak, radio features in galaxy clusters are categorised either as radio haloes or relics (Briiggen 
et al. 2012 ; Ferrari et al. 2008 ). Radio haloes are relatively symmetric, usually unpolarised, and 
cover the central parts of galaxy clusters, on scales of typically 1-2 Mpc. There are two main 
hypotheses to explain their origin: (re)acceleration of electrons up to relativistic energies due to the 
turbulence induced by cluster mergers (e.g. Brunetti et al. 2001 ), or collisions between thermal 
ions and relativistic protons in the ICM (Dennison 1980 ). Radio relics are typically elongated in the 
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Table 1 . Siizakii observation Logs 



Obs. ID 

Obs. Date 

R.A. 

Dec. 

Exposure 






(ks) 

W2 

ae808090010 

2013-07-02 

12 56 43.97 

-1-28 16 25.7 

12 

W3 

ae808091010 

2013-07-02 

12 56 15.34 

-1-28 20 14.6 

15 

Wl 

ae802084010 

2007-06-21 

12 57 22.27 

-1-28 08 25.1 

29 

W4 

ae806025010 

2011-12-09 

12 55 53.76 

-1-28 30 38.9 
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tangential direction and located near the clusters’ periphery. Their origin is thought to be connected 
to particle acceleration at shock fronts driven by the merger (e.g. Akamatsu & Kawahara 2013a , 
Ogrean et al. 2013a). 

Multi-wavelength studies of the thermal and non-thermal contents of clusters of galaxies us¬ 
ing X-ray, SZ, and radio data can therefore reveal important information regarding the mechanisms 
through which particles are accelerated (and the ICM is heated) during major cluster mergers. 

As one of the nearest, brightest, most massive merging systems, the Coma Cluster has been 
the target of a wealth of observations covering a wide range of the electromagnetic spectrum. The 
cluster is known to host both a giant radio halo (Coma C) and a radio relic towards the south-west 
(Willson 1970 ; Ballarati et al. 1981). X-ray observations have revealed the presence of a shock front 
associated with the radio relic (Akamatsu et al. 2013b ; Ogrean et al.2013b ). The giant radio halo 
exhibits a relatively sharp edge observed at 352 MHz with the Westerbork Synthesis Radio Telescope 
(WSRT, Brown & Rudnick 2011). 

Using SZ observations of the Coma Cluster, Planck Collaboration et al. 2013 indicated the 
existence of pressure jumps toward the south-west and south-east of the Coma Cluster centre, whose 
positions coincide with the edge of the WSRT giant radio halo. These pressure jumps are also asso¬ 
ciated with jumps in the X-ray temperature, especially along the western azimuth, and are therefore 
interpreted as weak shocks. The Mach number calculated from the pressure discontinuity along the 
western direction reported by Planck Collaboration et al. 2013 is Adwest = 2.03lo;o4. Using Suzaku 
observations, Simionescu et al. 2013 also report the detection of a temperature jump likely associated 
with a weak shock along the western side of the Coma Cluster, however located at a different radius 
and azimuth than the shock reported by Planck Collaboration et al. 2013 . 

In this paper, we present the results from two new Suzaku observations located along the west¬ 
ern edge of the Coma radio halo, which were obtained in addition to the data presented in Simionescu 
et al. 2013 . We have assumed a A cold dark matter cosmology with = 0.27, Ua = 0.73 and 
Ho = 70 km s“^Mpc“^. The Coma Cluster redshift is ^ = 0.0231 (Struble & Rood 1999 ), and one 
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arcmin corresponds to 28 kpc. We define the virial radius as r 2 oo> the radius within whieh the mean 
enelosed density is 200 times the eritieal density of the Universe at the redshift of the eluster. For 
the Coma Cluster, r 2 oo corresponds to 2 Mpe (70 aremin), with eonsistent values derived both from 
the Yx = Mgas X kT parameter (Planek Collaboration et al. 2013 ) and from weak leasing analysis 
(Okabe et al. 2014 ). 

2 Observations and Data Reduction 

We analyzed four Suzaku pointings loeated towards the west (W) from the Coma Cluster eentre, with 
a total exposure time of 73 ks. We denote these data sets as Wl, W2, W3 and W4, with inereasing 
index indieating a larger distanee of the aim-point with respeet to the eentre of the Coma Cluster. 
Wl and W4 were part of the mosaie analysed by Simioneseu et al. 2013 , while W2 and W3 are 
new observations performed by Suzaku in 2013 July (AO-8). The observation logs are summarized in 
Table 1. 

We redueed the data from the X-ray Imaging Speetrometers (XIS) 0, 1 and 3 using HEAsoft 
tools (version 6.13) and the ealibration database (CALDB) 2015-01-05. In addition to the standard 
sereening using the eriteria reeommended by the XIS team, we seleeted only observing periods with 
the geomagnetie eut-off rigidity COR > 6 GV. For the XISl data taken after 2011 June 1st, when the 
eharge injeetion level was inereased, we exeluded two adjaeent eolumns on eaeh side of the eharge 
injeeted eolumns. We also exeluded the eolumns in Segment A of XISO that were affeeted by an 
anomaly in 2009 June. 

Redistribution Matrix Files (RMF) and Aneillary Response Files (ARF) were ereating using 
the FTOOLS tasks xisrmf gen and xissimarf gen. 

3 Background Analysis 

Precise background subtraction is important for the spectral analysis of cluster outskirts, where the 
surface brightness is low. We, thus, carefully estimated both the non-X-ray background and cosmic 
X-ray backgrounds as described below. 

3.1 Non-X-ray Background 

The non-X-ray background (NXB) spectrum corresponding to each detector and each spectral extrac¬ 
tion region was created from night Earth observations using the task xisnxbgen. Tawa et al. 2008 
demonstrated that xisnxbgen reproduces the NXB spectrum within a 5 % accuracy. 

In order to further check the reliability of the NXB spectra for the newest observations (w2 
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Fig. 1. Surface brightness map of the four Suzaku pointings analysed in this work, using the 0.7-7.0 keV energy band. The image has been corrected for 
instrumental background and vignetting. White circles indicate the point sources excluded from the analysis. 352 MHz WSRT radio contours from Brown & 
Rudnick 2011 are shown in red. The contours start at 4mjy/beam and increase in intervals of 3 mjy/beam. The magenta circle marks a distance of 40 
arcmin from the Coma Cluster centre, while the yellow circle denotes a distance of 40 arcmin from the "analysis centre” chosen to match the radio halo edge. 


and w3), we employed the method suggested by Seta et al. 2013 . Briefly, for the front-illuminated 
(FI) devices (XIS 0 and 3), we measured the strength of the Ni I K a instrumental line at 7.47 keV 
in the NXB and the observed (non-background corrected) spectra by fitting these with a power-law 
plus Gaussian model in a narrow energy band centred around this emission line. We find that the 
normalisations of the Gaussian model agree between the observed and the corresponding NXB spec¬ 
tra, with ten out of twelve data sets (for the extraction annuli described in Section 4.2) showing an 
agreement within the la confidence interval, and two data sets showing an agreement between l-2cr, 
in accordance with the expected statistical dispersion. For the back-illuminated device (XISl), which 
does not show instrumental emission lines at high energies, we compared the count rate between the 
observed and NXB spectra in the 8-10 keV range and also find a good agreement given the expected 
statistical dispersion. We concluded that the NXB files created by xisnxbgen are a good description 
of the instrumental background corresponding to our observations. 

3.2 The X-ray backgrounds 

The X-ray background components are the Cosmic X-ray Background (CXB) due to the sum of 
distant, unresolved point sources, the Galactic Halo emission (GH) and the Local Hot Bubble (LHB). 
We adopt the same background model parameters determined by Simionescu et al. 2013 , using 
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regions in the Suzaku mosaie observations of the Coma Cluster that were loeated just beyond its 
virial radius. The effeet of Galaetie absorption was aeeounted for using the phabs model with the 
eolumn density iVn fixed at the average value eomputed from the Leiden-Argentine-Bonn radio HI 
survey (Kalberla et al. 2005 ) for the eoordinates eorresponding to eaeh Suzaku field of view. 

The CXB was modelled as a single power-law modified by the Galaetie absorption. The power 
law index and normalization reported by Simioneseu et al. 2013 are Tcxb = 1.50 and NormcxB = 
1.15 X 10“^ keV“^cm“^s“^ per 20^7r aremin^ at 1 keV, respeetively. The GH was modelled as an 
absorbed thermal plasma with temperature /cTgh = 0.25 keV, normalisation Normcu = 2.9 x 10“^, 
and 0.3 Solar metallieity, while the LHB model eonsisted of an unabsorbed thermal plasma with 
kTi^uB = 0.10 keV, Normi^uB = 1.6 x 10“^, and Solar metallieity assuming the abundanee table of 
Grevesse & Sauval 1998 . The normalizations of the thermal eomponents eorrespond to / neUndV x 
(10“^'^/(47r[i9A(l + ^)]^)) cm“® per 20^7r aremin^ , where Da is angular diameter distanee, rzg and 
rill are eleetron and hydrogen densities, respeetively, and 2 ; is the redshift. 

4 Analysis and Results 

To investigate the shoek along the western azimuth, we earned out temperature and surfaee brightness 
analyses, eonsidering the energy range between 0.7 — 7 keV. 

4.1 X-ray images 

To obtain a mosaie image of the region of interest, we eombined the images from all XIS deteetors 
and all pointings listed in Table 1. We exeluded anomalously bright pixels as well as the regions 
loeated within 30 aresee from the deteetor edges. Images of the instrumental baekground (NXB) 
sealed to the eorresponding exposure time for each detector and each observation were subtracted 
from the raw count maps. To take into account the vignetting effect, we created flat-field maps using 
the task xissim, which performs a ray-tracing Monte Carlo simulation of the X-ray telescope (XRT) 
and XIS^ The input photon list for xissim was that representing a spatially uniform emission with 
a spectrum corresponding to the best-fit average model of the cluster emission and sky background 
obtained from the central 6.5' of each respective pointing. Using this model spectrum, the ray trac¬ 
ing then accounts appropriately for the energy dependence of the vignetting effect. The vignetting 
corrected image was created by dividing the background subtracted XIS image by the flat-field map. 
The obtained image is shown in Figure 1. In the same figure, we indicate the 352 MHz WSRT radio 
contours (Brown & Rudnick 2011 ) in red. 

‘ see https://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/expomap.html 
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For the temperature and surface brightness profiles presented in the following sections, we 
have considered annuli with two different centres: one is the cluster centre, which we adopt as 
(a, 5) = (12h59m42.44s,+27°56'45.53"), and the other is referred to as the “analysis centre” lo¬ 
cated at (a, (5) = (13h00ml8.043s,-|-28°06'46.07"). The analysis centre was defined in such a way 
that a circle of radius 40 arcmin centred on this point traces most closely the western edge of the 352 
MHz WSRT radio contours (see yellow dotted curve in Figure 1). 

4.2 Spectral Analysis 

To derive the temperature profile, we extracted and modelled spectra from six annuli with radii 26- 
30, 30-35, 35-40, 40-45, 45-52, 52-60 arcmin from the cluster centre, and 30-35, 35-40, 40-45, 
45-50, 50-55, 55-62 arcmin from the “analysis centre”. Visually identified point sources or compact 
substructures marked by white circles in Figure 1 were excluded from the spectral extraction regions. 

The spectral modelling was performed with the XSPEC spectral fitting package (version 
12.8.2, Arnaud 1996 ), employing the modified C-statistic estimator. Unless otherwise noted, sta¬ 
tistical errors are reported at the AC = 1 confidence level. The spectra were binned to a minimum 
of one count per channel and fitted in the 0.7-7.0 keV band. As the spectral model, we employed an 
absorbed thermal plasma in collisional ionization equilibrium using the apec code (Smith et al. 2001 
), in addition to the X-ray background model described in Section 3.2. 

The free parameters in the spectral model are the normalization and temperature of the apec 
component for each annulus. The apec metal abundance was treated as a free parameter but coupled 
to have the same value for all annuli. We obtain a best-fit value of Z = 0.20 ± 0.07 Solar in the units 
of Grevesse & Sauval 1998 . Within the la statistical confidence interval, this value is consistent with 
that reported by Werner et al. 2013 for the outskirts of the Perseus Cluster. 

4.3 Temperature profile 

The values for the best-fit parameters obtained from the spectral fitting are summarised in Table 2. 
Figures 2 and 3 show the radial profiles of the temperature obtained from this analysis for the two 
assumed sets of annuli. 

The temperature profile obtained from annuli centred on the cluster centre shows a marked 
temperature jump between the annuli spanning 35 — 40 arcmin and 40 — 45 arcmin (Figure 2). This 
indicates the presence of a shock front. 

From the Rankine-Hugoniot shock jump condition (see Landau & Lifshitz 1959 ; Markevitch 
& Vikhlinin 2007 ), the Mach number can be estimated from the temperature discontinuity using the 
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Table 2. Temperatures, normalizations, and column densities of the Galactic absorption for each annulus. 


The 


unit of the spectrum 


normalization Is defined as f neriHdV x (10“^''/(47r[DA(l + 2 )]^)) cm“® per 20^7r arcmln^. 


Radius (aremin) 

kT (keV) 

Norm (X10 

(Vh (lO^^* cm-2) 

centering on the cluster centre 

26-30 

Q IQ+0-71 

4.40 4=0.07 

8.73 

30-35 

8.672“;^i 

2.47 4=0.04 

8.73 

35-40 

lo.ertlir 

1.20 =h 0.03 

8.61 

40-45 

^.O^_0 Q2 

0.69 4=0.04 

8.48 

45-52 

r on+0.91 

o.ot;_Q Y3 

0.67 4=0.03 

8.48 

52-60 

O.OU_Q gj 

0.262^1 

8.48 

centering on the analysis centre 

30-35 

8.592“:«° 

4.614=0.07 

8.73 

35-40 

9.402“/j 

2.814=0.04 

8.73 

40-45 

8-9421:“ 

1.28 4:0.03 

8.61 

45-50 

r qq+114 
tJ.OO_o.gi 

0.60 4=0.03 

8.48 

50-55 

r. 97+0.91 
( -0.65 

0.58 4=0.02 

8.48 

55-62 

+.DO_o.65 

0.42 4=0.02 

8.48 


following equations: 

/I 1 

where 7 = 5/3 is the adiabatie index, p = pi/po is the eleetron density jump, t = Ti/Tq is the gas 
temperature jump, and C = (7 + l)/(7 — !)■ The indiees 0 and 1 indieate the values of the respeetive 
quantities before and after the shoek passage, respeetively. Based on the observed temperature jump 
and taking into aeeount the propagation of errors, we estimated the Maeh number as Ad = 2.3 ± 0.4. 

However, the Maeh number determined in this way is only valid under the assumption that the 
entirety of the temperature differenee between the annuli at 35 — 40 aremin and 40 — 45 aremin is due 
to shook heating or - otherwise put - that no temperature gradient between these regions was present 
before the shook passage. In reality, the temperature is expeeted to deeline with radius in the oluster 
outskirts, even in the absenoe of a shook. From iV-body simulations. Burns et al. 2010 prediot that, 
for relaxed elusters, the temperature profile as a funotion of radius oan be desoribed by the following 
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Fig. 2. Suzaku temperature profile for the Coma Cluster. Black points were obtained in this work. The other measurements, shown with dashed lines, are 
reproduced from Simionescu et al. 2013 . Green, red and blue represent the west, north west and east directions, respectively. The gray curve indicates the 
universal profile expected from numerical simulations of relaxed clusters (Burns et al. 2010 ) scaled to the average temperature of this system. {R 200 = 70 
arcmin.) 


equation: 


= A 


-avg 


1 + B 


/ j. 


200 


1-/3 


(3) 


where Tavg is the average X-ray weighted temperature between 0.2 and 1.0 r 2 oo- Bums et al. 2010 
report the best-fit values of the other parameters as 4 = 1.74±0.03, i? = 0.64±0.10 and /3 = 3.2 ±0.10. 
Using the measurements along all azimuths of the Coma Cluster covered by Suzaku (this work and 
Simionescu et al. 2013 ), we obtain Tavg = 8.43 ± 0.05 keV, resulting in the predicted temperature 
profile shown as a grey band in Figure 2. The measured temperature jump at the 40 arcmin shock is 
clearly higher than the temperature gradient expected for a relaxed cluster. The pre-shock temperature 
in the 35 — 40 arcmin annulus predicted from the universal temperature profile of Bums et al. 2010 
is Tuniv. = 5.71 ± 1.0 keV, which would imply a Mach number Af = 1.8 ± 0.4 - smaller than, but 
roughly consistent with the value derived without accounting for the underlying temperature gradient 
in the absence of a shock. 

We also evaluated the temperature profile from annuli centred on the “analysis centre” (see 
Figure 3). In this case, the temperature jump appears between the annuli at 40-45 and 45-50 arcmin, 
and the corresponding Mach number is estimated as A4 = 1.7±0.3. Note that, in this case, due to 
the special geometry, we cannot apply the same correction for the underlying temperature gradient 
before the passage of the shock as discussed in the previous paragraph and our value is therefore an 
upper limit for the measured Mach number. 
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Fig. 3. Temperature profile centred on the analysis centre. A discontinuity can be confirmed around 45 arcmin. The green line indicates the location of the 
radio halo edge. 

4.4 Systematic uncertainties on the temperature estimation 

In order to account for the systematic uncertainties associated with the X-ray background model, we 
divided the area beyond the cluster’s r 2 oo from which the average model parameters were estimated by 
Simionescu et al. 2013 into 16 different spectral extraction regions, and fitted each region separately 
with the normalizations of the CXB power-law and GH thermal emission allowed to vary. All other 
parameters of the X-ray background model were kept fixed to the values described in Section 3.2. 
Each such background region roughly corresponded to a half of the Suzaku XIS field of view. 

We then calculated the mean and standard deviations of the CXB and GH normalisations from 
these 16 different regions and refitted the radial temperature profile by setting each parameter in turn 
to its mean—Icr and mean-flcr values. In practice this amounts to a ±9% variation of the CXB and 
±30% variation of the GH fluxes. Fig. 4 shows the results. We conclude that the effect of systematic 
uncertainties is small. This is not unsurprising, since the shock is located at radii of only about 40 
arcmin, less than two-thirds of the Coma Cluster’s r 2 oo- 

4.5 Surface brightness profile 

The presence of a shock feature in the ICM should cause not only a temperature jump, but also a 
discontinuity in the gas density and therefore the X-ray surface brightness. To search for such a 
density jump associated with the shock, we extracted surface brightness profiles from the XIS mosaic 
shown in Figure 1, using the energy band between 0.7 and 7.0 keV. 

We then used the PROFFIT code vl.l (Eckert et al. 2011) to fit the surface brightness profiles 
with a projected broken power-law model, assuming spherical symmetry with respect to the centre of 
the annuli used to obtain the radial profile: 
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Fig. 4. The expected effects of spatial variations of the CXB and GH normalisations on the measured radial temperature profile. The blue and red lines 
indicate temperature changes by varying the CXB and the GH normalizations within their la- confidence intervals, respectively. 


— C-Tlnorm f ^ , T < ?"break (4) 

X ^break X 

^0 ^norm ( ) , T" > ^"break- (5) 

X ^break X 

In the above equations, n is the eleetron number density, C is the density eompression faetor, r is 
the distanee from the centre of the sector, rbreak is the break radius and a is the power-law index. 
Subscripts 1 and 0 indicate the post-shock and pre-shock properties, respectively. 

To obtain the fitting results, we adopted the Cash statistics, which is appropriate for the low- 
count regime. We include the sky background as an additional constant surface brightness component 
in our model. The value of this constant was estimated from the same Suzaku pointings that were 
used to determine the X-ray background model used for spectral analysis (these background fields are 
located beyond the 70 arcmin, the cluster’s estimated r 2 oo)- We include the effects of Suzaku’s point 
spread function (PSF) by convolving the surface brightness model with a Gaussian with a full width 
at half maximum of 1.9 arcmin. 

We first fitted the radial surface brightness profiles centred on the cluster centre, restricting 
our fit to the radial range between 26 and 45 arcmin. The result is shown in Figure 5. The best 
fit parameters are ao = 0.76lo;og, ai = 2.59 ± 0.13, rbreak = SS.dlgJo arcmin and C = 1.31 ± 0.06. 
Accordingly, the corresponding Mach number from the Rankine-Hugoniot condition for the density 
jump is Ad = 1.21 ±0.04. 

For fitting the surface brightness profiles centred on the “analysis centre”, we used the radial 
range between 29 and 50 arcmin. The density discontinuity is stronger and more easily seen, and 
the computed best fit parameters are ao = O.dllo;!]®, ai = 2.59llo^l, ''"break = 40.5^0 2 ? arcmin and 
C = 1.65lo 09 , resulting in the Mach number Ad = 1.45 ±0.08. As can be seen in the right-hand panel 
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Fig. 5. Surface brightness profiles of the Coma Cluster and the best-fit projected broken power-law models centred on the cluster centre (left) and the “analysis 
centre” chosen to match the curvature of the outer 352 MHz WSRT radio contours (right). In both panels, the red lines indicate the broken power-law model 
and the blue lines indicate the instrumental background, NXB. The green vertical shaded region in the left panel indicates the radial range of the radio halo 
edge. In the right panel, a dark-green line indicates the surface brightness profile obtained from the 352 MHz map. 


of Figure 5, the location of the edge in the X-ray surface brightness profile computed with respect to 
the “analysis centre” corresponds well to the location of the radio halo’s edge. 

We have also searched for spatial variations of the Mach number by calculating surface bright¬ 
ness profiles centred on the “analysis centre” in two separate wedges with equal opening angles. We 
find no evidence for significant differences in the shock strength, with rbreak = 40.8 ± 0.4 arcmin and 
A4 = 1.47 ± 0.15 for the northern half and rbreak = 40.4 ± 0.5 arcmin and A4 = 1.32 ± 0.10 for the 
southern half of the shock covered by Suzaku. 

5 Discussion and Conciusions 

We have analysed data from four Suzaku pointings located along the western side of the Coma Cluster. 
Both the temperature and surface brightness distributions reveal the presence of a shock front likely 
associated with the edge of the Coma radio halo. In this section, we investigate the strength and 
location of the shock, in comparison with previous results estimated from X-ray, SZ, and radio ob¬ 
servations, and discuss particle acceleration at X-ray shocks in relation to the origin of the radio 
emission. 

5.1 Shock strength, location, and geometry 

Previous X-ray studies have already revealed the presence of temperature structure in the Coma 
Cluster. ASCA observations (Honda et al. 1996 ; Watanabe et al. 1999 ) showed a temperature 
asymmetry, with the western side of the cluster being hotter than the eastern side. This trend was later 
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confirmed by XMM-Newton and Suzaku (e.g. Amaud et al. 2001 ; Simioneseu et al. 2013 ). This 
temperature strueture is thought to be assoeiated with a past merger (Honda et al. 1996 ; Amaud et 
al. 2001 ). 

Using SZ observations, Planek Collaboration et al. 2013 report the presenee of two shook 
fronts looated about 30 aromin towards the west and south-east from the eluster oentre, with Maoh 
numbers A4w = 2.03lo;o4 and A4se = 2.05lo;o2 derived from the Rankine-Hugoniot pressure jump 
oonditions. Both of the shooks reported by Planek Collaboration et al. 2013 align well with the 
sharp gradients in the 352 MHz WSRT radio brightness distribution, and are assoeiated with X-ray 
surfaoe brightness edges and temperature jumps deteoted in XMM-Newton data. Thus, shook heating 
is likely at least partly responsible for the elevated gas temperatures on the western side of the eluster 
reported by previous studies. 

In this work, we also report the presenee of a shook front along the western side of the Coma 
Cluster, although looated further north than the western shook reported by Planek Collaboration et al. 
2013 ; the area oorresponding to the shook seen by Planek is not oovered by Suzaku observations. For 
different ohoioes of extraotion regions and different assumptions about the underlying temperature 
profile in the absenee of shook heating, the temperature profiles measured with Suzaku imply Maoh 
numbers ranging from A4 = 1.7 ± 0.3 to Ad = 2.3 ± 0.4, broadly oonsistent with the Maoh number 
derived from SZ observations for what oan be assumed to be a different azimuth along the same shook 
front that traoes the entire western edge of the giant radio halo. 

Shook signatures are seen not only in the temperature but also in the X-ray surfaoe brightness 
distribution along this azimuth. When extraeting radial profiles in annuli roughly parallel to the steep 
radio edge, we obtain oonsistent Maoh numbers from the inferred temperature and surfaoe brightness 
jumps, with Ad = 1.7 ± 0.3 and A4 = 1.45 ± 0.08, respeotively. However, when ehoosing annuli een- 
tred on the peak X-ray emission of the Coma Cluster, the shook strengths derived from the tempera¬ 
ture and surfaoe brightness beeome ineonsistent, with the Maoh number inferred from the temperature 
jump. Ad = 2.3 ± 0.4, being larger than that determined from the density jump. Ad = 1.21 ± 0.04. 
Even when aeeounting for the expeeted underlying radial temperature gradient between the pre- and 
post-shook regions in a typioal relaxed eluster, the Maoh number inferred from the temperature jump 
remains high, at Ad = 1.8 ± 0.4. 

It is generally expeeted that a mis-alignment of the shook front eompared to the orientation 
of the speotral extraotion regions would lead to an underestimation of both the true temperature as 
well as density jumps. Therefore, the observed deerease of the Maoh number estimated from the 
surfaoe brightness profile, from Ad = 1.45 ± 0.08 to 1.21 ±0.04 for a less optimised geometry, is 
easily understood; however, assuming the shook has a spherioally symmetrioal shape, it is diffieult to 
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explain why the shoek strength inferred from the temperature jump does not show a eommensurate 
deerement. 

Perhaps even more surprisingly, the loeations of the diseontinuities determined from the X-ray 
temperature and surface brightness differ. When considering the profiles centred around the “analysis 
centre”, we obtain rbreak = 40.5loJi arcmin from the surface brightness profile fitting; this is in 
good agreement with the location of the radio halo edge (the “analysis centre” was chosen so that 
a circle of radius 40 arcmin centred on this point traces most closely the western 352 MHz WSRT 
radio contours). On the other hand, the temperature jump is seen at a radius of 45 arcmin; a quick 
inspection of Figure 3 reveals that, for a shock located at a radius of 40.5 arcmin, we expect to detect a 
temperature jump between the second and third annuli (35-40 and 40-45 arcmin), rather than between 
the third and fourth annuli as currently seen. A similar conclusion is reached for profiles centred 
around the cluster centre. In this case, the surface brightness jump is seen at 33.4lg gg arcmin, in good 
agreement with the location of the radio edge contours spanning 34-35 arcmin with respect to this 
choice of centre, while the temperature jump is seen at 40 arcmin. 

Thus, the locations of the temperature and density jumps appear to differ by 4.5-6.6 arcmin 
(125-185 kpc), with the temperature jump located further out while the surface brightness jump agrees 
well with the position of the edge of the radio halo. This difference is significantly larger than the 
Suzaku PSF, and is seen regardless of the choice of the centre used to extract the profiles. If the 
locations of the temperature and density jumps are better aligned for the southwestern half of the 
region characterised by a steep radio gradient, this would explain why the northwestern shock has not 
been previously identified either by Planck Collaboration et al. 2013 (the pressure jump is smeared 
because the product of temperature and density now has discontinuities at multiple radii), or by Brown 
& Rudnick 2011 who showed that XMM-Newton spectra from smaller extraction regions on either 
side of the radio edge only suggest a significant temperature jump towards the southwest but not 
towards the northwest. 

A possible explanation for the different shock radii inferred from the temperature and surface 
brightness profiles is that the shocked gas that covers the largest line of sight depth (and has con¬ 
tributed to accelerating the most radio-emitting particles) has been heated by a weaker shock and 
lies closer to the centre, while the temperature jump is driven by gas affected by a stronger shock 
at slightly larger radii, but where the shock has a smaller line of sight depth; while this might still 
steepen the power-law fitted to the surface brightness distribution, the line-of-sight effect may con¬ 
ceal an actual discontinuity corresponding to the location of the temperature jump. We note that 
numerical simulations do indeed suggest both complex shapes as well as complex distributions of 
the Mach number along the shock fronts believed to be responsible for accelerating particles to the 
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relativistic energies required to produce radio emission in radio relics (Skillman et al. 2013 ), so it 
is plausible that this may also apply for the Coma radio halo. An alternative, perhaps more exotic 
scenario, is that in which a fraction of the particles accelerated by the shock may diffuse ahead of the 
shock and cool out of the radio band by mixing with and heating the thermal ICM. 

5.2 Particle acceleration at X-ray shocks and the origin of the radio emission 

The two main hypotheses typically invoked in the past in order to explain the synchrotron emission 
in radio haloes are typically (re)acceleration of electrons up to relativistic energies due to the turbu¬ 
lence induced by cluster mergers (e.g. Brunetti et al. 2001 ), or collisions between thermal ions and 
relativistic protons in the ICM (Dennison 1980 ). Concomitantly, in a growing body of publications, 
the edges of radio haloes are found to be coincident with shock fronts (e.g. Macario et al. 2011 , 
Markevitch et al. 2005 ; Shimwell et al. 2014 , Planck Collaboration et al. 2013 , and this work). This 
suggests that shocks are an additional process contributing to the acceleration of particles in radio 
haloes. The clusters where shocks are seen to be associated with the edges of radio haloes therefore 
present a unique angle to study a superposition of particle acceleration mechanisms acting simultane¬ 
ously. The Coma Cluster is the nearest, brightest, and to date one of the best studied among this type 
of systems. 

One of the main theoretical challenges to be addressed in this case is the exact mechanism 
through which very weak shocks (Af ~ 2 in Coma and typically < 3 in other clusters cited above) 
are in fact able to boost the radio surface brightness. For such low Mach numbers, the diffusive 
shock acceleration (DSA) mechanism alone is thought to be inefficient at accelerating particles. The 
preferred interpretation is that, instead of directly accelerating the thermal electrons, the shocks re¬ 
accelerate a population of already mildly relativistic electrons present in the cluster (e.g. Markevitch 
et al. 2005 , Kang & Ryu 2011, Pinzke et al. 2013 ). 

The morphology of the Coma radio halo supports this interpretation; the observed steepening, 
or “edge”, in the radio surface brightness profiles reported by Brown & Rudnick 2011 , is sharpest 
along the western side of the cluster, and is not apparent along most other directions. Suzaku demon¬ 
strates for the first time that a shock front extends not only along the southern but also along the 
northern half of this radio brightness edge. Therefore, we can infer that the difference in radio bright¬ 
ness between the western side of the cluster and other azimuths is entirely attributed to reacceleration 
due to the shock passage. 

Recent particle-in-cell simulations indicate that a different mechanism termed shock drift ac¬ 
celeration may result in a more efficient acceleration at low-Mach number shocks compared to stan- 
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dard DSA (e.g. Caprioli & Spitkovski 2014 ; Guo et al. 2014a ; Guo et al. 2014b ), further boosting 
the radio surfaee brightness at the western Coma shook. 
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Fig. 6. Observed XIS spectra for each annulus centered on the cluster center. (1) 26-30 arcmin (2) 30-35 arcmin (3) 35-40 arcmin (4) 40-45 arcmin (5) 
45-52 arcmin (6) 52-60 arcmin. We model the spectra with a thermal component describing the ICM of the Coma Cluster (Section 4.2) in addition to a model 
describing the X-ray backgrounds (Section 3.2). The solid lines indicate the best fit model. We show separately the data from all three detectors and different 
overlapping pointings within each annulus: XISO (black and blue), XIS1 (red and cyan) and XIS3 (green and magenta). 
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